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Charge and spin current correlations are analyzed in a source of spin-entangled electrons built 
from a superconductor and two quantum dots in parallel. In addition to the ideal (crossed Andreev) 
channel, parasitic channels (direct Andreev and cotunneling) and spin flip processes are fully de- 
scribed in a density matrix framework. The way they reduce both the efficiency and the fidelity of 
the entangler is quantitatively described by analyzing the zero- frequency noise correlations of charge 
current as well as spin current in the two output branches. Spin current noise is characterized by a 
spin Fano factor, equal to (total current noise) and —I (crossed correlations) for an ideal entan- 
gler. The violation of the Bell inequalities, as a test of non-locality (entanglement) of split pairs, 
is formulated in terms of the correlations of electron charge and spin numbers counted in a specific 
time window r. The efficiency of the test is analyzed, comparing r to the various time scales in the 
entangler operation. 
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I. INTRODUCTION 

The production and the analysis of entangled states 
in condensed matter devices have recently emerged as 
a mainstream in nanoelectronics. Indeed, entangle- 
ment between electrons, besides checking fundamental 
quantum properties such as non-locality, could be used 
for building logical gates and quantum communication 
devices*. One may choose to consider nanoscopic de- 
vices where electrons behave essentially as free particles 2 , 
in a way similar to photons in quantum optics. Alter- 
natively, exploiting electron interactions open new pos- 
sibilities. This allows for instance to use the electron 
spin as a qubit in quantum dots*, owing to spin relax- 
ation/coherence times ranging in semiconductors from 
fractions of /is (bulk*) to fractions of ms (Ti in quan- 
tum dots*). Spin entanglement 6 jL* and teleportation 9 
scenarios have been proposed within such a framework. 
A test of this spin entanglement was previously proposed 
with the help of a Bell diagnosis^* , for a normal metal- 
superconducting device described within the context of 
scattering theory, operating in the ideal crossed Andreev 
regime 11 . However, in any device several unwanted elec- 
tronic transitions may spoil entanglement and introduce 
decoherence. It is of key importance to check the device 
robustness against these parasitic effects and to show how 
to diagnose them through transport properties. Owing 
to the strategic role of entangled pairs in quantum in- 
formation, we focus here on the source of spin-entangled 
electrons proposed by Recher et al£ Following a de- 
tailed study of the average current due to the various 
processes*^, the present work addresses noise correlations 
as a diagnosis of the device operation. 

The entangler of Ref. is made from a superconductor 
(S), adjacent to two small quantum dots in parallel, each 
being connected to a normal lead featuring a quantum 
channel. The dots are assumed to filter electrons one 



by one in a single orbital level, and Cooper pairs emit- 
ted from the superconductor are split in the two dots 
- the so-called crossed Andreev channel (CA)****. Ide- 
ally, the constituent electrons then propagate in the two 
output leads as an entangled pair. Parasitic quantum 
processes may also occur in the entangler, they include 
elastic cotunneling (CT)**** between the two dots via 
S, and tunneling of a pair through the same lead by a 
Direct Andreev process (DA) (Fig. The main (CA) 
channel as well as DA were studied in Ref. 0. Using 
perturbation theory, the average current contributions 
to CA and DA were separately calculated. In the same 
spirit, Ref. used a beam splitter on the entangler's 
output, in order to detect the entangled singlets by their 
noise bunching properties^*. In reality, all these trans- 
port channels arc mixed together. Instead of analyzing 
them separately, the full electron flow should be studied 
sclf-consistently, and be eventually characterized by its 
density matrix. The present work is based on a quantum 
master equation scheme which treats all processes on an 
equal footing** and gives access to the full density ma- 
trix of emitted electrons. It is used to compute the charge 
and spin correlations in the current flowing through the 
two branches of the entangler. In particular, spin current 
correlations* 9 ", even taken at zero frequency, are shown to 
quantify the entangler's efficiency and fidelity. Next, the 
violation of Bell inequalities (BI) can be tested by com- 
puting the cross correlators of charge and spin particle 
numbers in a given time window. 



II. IDEAL REGIME AND THE OTHER 
TRANSPORT PROCESSES 

Transport is described by the processes connecting the 
different states - identified by the charge number and 
spin of the dots. Extra electrons n±, Ti2 define the charge 
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FIG. 1: Processes involved in the operation of the entan- 
gler. Pairs are emitted by the superconductor (gray) towards 
the two quantum dots, electrons are then collected by the 
leads L, R along predefined spin polarizations ±a, ±6 (4" 1 
panel). The order of elementary transitions involving virtual 
states is indicated : Crossed Andreev (CA) produce split pairs 
and triggers entanglement. CT (cotunneling), Direct Andreev 
(DA) as well as one-by one tunneling through the same dot 
spoil entanglement. The three left panels describe coherent 
transitions, while the two right ones describe incoherent tran- 
sitions. 



states (n\n2) in the dots. The electronic processes at 
work in the device include both coherent processes (CA, 
CT and DA) and incoherent processes (transitions be- 
tween the dots and the leads), see Fig. ^ In the ideal 
operation of this entangler, the Coulomb blockade pre- 
vents double occupancy in each dot, so that DA is for- 
bidden. Also, the superconducting gap is supposed to 
be very large, so that one-by-one tunneling (last panel in 
Fig. ^) is also excluded. In addition, CT is neglected in 
this ideal regime, owing to the energy difference between 
the electron states on the two dots. Then, starting from 
an empty orbital state (denoted (00)), crossed Andreev 
(CA) reflection couples to a singlet state shared between 
the two dots, (ll s ) with an amplitude jaT (ja <C 1 is 
the geometrical factorii^i^) . For a resonant CA transi- 
tion, the dot energy levels Ei satisfy E\ + E 2 = 0. The 
two electrons are collected into the leads (with chemical 
potentials hl,r < E\, E 2 ) by transitions to states (01) or 
(10), then (00), with rates Ti (i = 1,2). 

Let us now consider parasitic transport processes. If 
the intradot Coulomb charging energy U is not very 
strong, a coherent transition from (00) to a doubly oc- 
cupied dot state (20) or (02) can occur via a direct An- 
dreev (DA) process^. Notice that this occurs with an 
amplitude T larger than for CA, since it is not affected 
by the geometrical reduction factor 7^4. Electrons are 
successively collected by a single lead, with rates T ■ then 
Ti, eventually reaching the empty state (00). This trans- 
port channel implies dot energies 2Ei + U. In addition, 
DA can also start from an initial state (10) or (01), pro- 
ceeding through states (21) (respectively (12)) with en- 
ergies 2Et +E 2 + U (respectively Ei + 2E 2 + U) . Further 
collection into the leads either give states (20), (02), or 
singlet and triplet states (ll Sl t). The mixing of triplet 
states with the desired singlets may cause decoherence in 
the source operation. In another process, the two elec- 
trons from a Cooper pair can tunnel one by one towards 
the same leadi, involving a singly occupied virtual state 
which costs an energy at least equal to the superconduct- 
ing gap As (Fig. 2J). Contrary to the DA process, the 



dot is emptied before the quasiparticle in S is annihilated. 
The state (10) (or (01)) is reached by an incoherent pro- 
cess with a rate Ti ~ TiT?/A 2 s . Last, CT involves a 
coherent transition of an electron from one dot to the 
other via S, and mixes all the above processes^: it cou- 
ples state (01) to (10), but also (20) (or (02)) to (11), 
(21) to (12). CT has an amplitude 7 C T < T, of the 
same order of magnitude as that of CA. 

The density matrix equations involving populations 
(diagonal elements) and coherences (nondiagonal ele- 
ments) include these processes altogether in a consistent 
and non-perturbative manner. Here one assumes weak 
couplings r^ij < fcs8 of the dots to the output leads L, 
R (0 is the temperature). On the other hand, the leads 
are biased such as |/i_L,j?J > T, Tl^r, ensuring irreversible 
collection of the electrons in L, R. The density matrix 
equations have been derived and detailed in Refs. ll2TIT^ . 
where the average currents in each branch have been com- 
puted. From these Bloch-like equations, one can also cal- 
culate the conditional probabilities - for a transition to 
a given state, assuming a previous initial state - which 
enter the calculations of the noise correlators22i2ii22 . Pro- 
cesses involving states with at most two electrons on the 
double dots are represented on Fig. 6 of Ref. 0, yet 
three-electron states are also included in the calculations. 



III. CHARGE AND SPIN CURRENT 
CORRELATIONS 

Let us assume that the leads contain separate channels 
perfectly filtering spin currents polarized along two pre- 
defined and opposite directions. An understanding of the 
charge and spin correlations can be obtained by consid- 
ering the spin-resolved zero-frequency noised : I a being 
the current carried by electrons with spin a = ±a (along 
a given spin direction a), we define 



(0) 



dt({AI?(t),AI?(0)}) , 



(1) 



where i,j = {L,R} and AI[(t) is the deviation from 
the average current component (If(t)) = If . We assume 
that the lead resistances and couplings to the dots are not 
spin-dependent. Because the superconductor is emitting 
singlet pairs, it results that the average current in each 
lead is not spin-polarized, e.g. I? — h/2. By definition, 
the charge and spin current noise correlation read: 
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The latter simply expresses the time correlations of the 
spin current I^ pin = I? — lf a . The average total charge 
current passing through the entangler can be written as 

I = Il+ Ir, (4) 
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where the currents in leads L and R can be separated in 
two components, as 

h = + Ilr , (5) 

Ir = 2I RR + I LR . (6) 

Here III, Irr and Ilr respectively count pairs pass- 
ing in unit time through L, or R, or through (L, R) as 
split (entangled) pairs. These current components can 
be written as: 

Ilr = \p S I , (7) 

III = \vlI , (8) 

Irr = t;Pr! , (9) 

with probabilities such that pl + pr + Ps — !■ 

It is crucial to notice here that all the electrons pass- 
ing through the entangler enter the two-dot system by 
pairs, emitted in an Andreev process (CA, DA or onc- 
by-one tunneling). Therefore, the spin current fluctu- 
ations due to electrons of a given pair are correlated, 
since they track the passage of singlets. On the contrary, 
electrons emitted within different pairs display no spin 
correlation whatsoever. It results that, while calculat- 
ing the correlations of spin currents at different times, 
the contributions of separate pairs drop out, so that spin 
noise in the leads uniquely tracks the spin correlations in- 
side a pair. These correlations should be non-ideal since 
they track what remains from the emitted singlets after 
their passage through the two-dot system. This property 
of probing single pairs strongly contrasts with the charge 
current fluctuations which - in a sequential tunneling pro- 
cess - correlate successive pairs due to the Pauli principle 
and the Coulomb interaction in the dots. Therefore the 
spin current noise is very well suited to study the spin 
correlations inside pairs with zero-frequency noise. To 
illustrate this, it is convenient to write the various spin 
noise correlations as follows : 

S? L = 2eI L - 4eI LL (l - f L ) = 2eI[p L f L + \ P s] , 

S s g R = 2eI R - 4eI RR (l - f R ) = 2eI\p R f R + \p s ] , 

Slr = S S rl = -2eWl - fs) = -elpsQ. - fs) ■ 

(10) 

Let us comment the various terms in Eq. (10). The first 
term in S^f corresponds to the autocorrelation of elec- 
tron wave packetsi2i21i ; and takes a Poisson value. The 
other terms are negative and reflect the spin anticorre- 
lation inside singlets. The Fano-like reduction factors 
fi, fs quantify this correlation. For instance, blocking 
one lead (Tr ; = for instance), all pairs pass through L 
(pl = I ) in a sequential way. In this case I = II = 21 ll- 
Moreover, since the dot filters electrons one by one, one 
easily checks that f L — 0, thus the spin noise through 
L is zero. This "spin Fano factor" equal to zero is a 
property of a S/N junction, characteristic of an unspoiled 



Andreev process^ : at a superconductor/normal (S/N) 
metal junction, zero spin noise is the fingerprint of charge 
carriers being paired in a singlet state. On the other 
hand, if the entangler operates ideally (only CA), then 
ps = 1) fs = an d — el. This is somewhat sim- 
ilar to the behaviour of a quantum dot in the sequen- 
tial regime, with a single orbital state involved in trans- 
port and strong Coulomb charging energy, yielding a spin 
Fano factor equal to onei^. 

On the other hand, the crossed spin correlation is 
S L P R — —el. This value corresponds to an effective spin 
Fano factor —1, and serves as a reference value for an 
ideal operation of the device. These Poisson-like Fano 
factors, confirmed - see below - by the solution of the 
master equations, are valid for any value of the couplings 
of the dots to the leads. This is in contrast with the 
"charge" Fano factors in resonant quantum dot devices 
which reach the value one only in the very asymmet- 
ric limit, where correlations between successive charges 
crossing the dots become negligible. 

Summing all contributions in Eq. (10), the spin noise 
of the total current is 

SZ = S l P l + S R P R + S? R + S% L = 2eIF s , (11) 
with the total spin Fano factor 

Fs = Vhfh + PR.fR + psfs ■ (12) 

Notice that in our definition the spin Fano factors mea- 
sure the zero-frequency spin current noise in units of 
the average charge current (not the average spin current 
which is zero). 

In the general case, when all (CA, DA and CT) pro- 
cesses come into play, pairs are both distributed (e. g. 
they exit through LL, RR and LR) and " mixed" together 
(electrons of a pair may intercalate between electrons of 
another) : for instance, due to states with three elec- 
trons in the double dot, two electrons of a split pair L, R 
can be separated by a time interval during which one 
or more pairs pass through a single dot (as illustrated 
by the sequence (00), (11), (10), (12), (11), (10), (00)). 
This variable delay between two spin-correlated electrons 
is responsible for the factors fL,R,s > which measure 
the degree of mixing. 

This preliminary analysis suggests that parasitic pro- 
cesses have two effects. First, they reduce the proba- 
bility of split (entangled) pairs over non-entangled ones, 
thus p s defines the entangler efficiency. And secondly, 
they reduce the spin correlation of a split pair, 1 — f s 
defining the entangler fidelity with respect to the singlet 
state. Notice that both efficiency and fidelity are defined 
here with respect to a given measurement probe, e. g. 
zero-frequency noise, or alternatively time-resolved cor- 
relation in the present work. 
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IV. QUANTITATIVE ANALYSIS 

We now study these trends quantitatively. Charge and 
spin current noise correlations are calculated from the 
entangler quantum master equations, derived in ref. 0, 
where four-electron states in the dot pair are neglected 
for simplicity. We successively focus on the correlations 
of the total current II + Ir, of the separate currents II, 
I R , and on the crossed correlations between L and R. 

T/T U/T 



1 2 3 4 10 20 30 40 50 




FIG. 2: Spin current cross-correlations, assuming CA to be 
resonant: E\ + E2 = 0, and Tl = Tr. As is the largest 
energy scale, thus one-by-one tunneling is negligible. T = 1 
and j a ~ 0.2 is assumed. In a) and b), only CA and CT 
are taken into account, assuming U —* 00. Charge and spin 
cross-correlations are plotted for < Y < 2 and several values 
of E = \Ei - E 2 \ = {0,0.2,0.5,1,5} (bottom to top). The 
dotted line features the ideal case, c) and d) correspond to 
< U < 50 and E = 5 (CA and DA, but no CT), for Y l ,r = 
{0.01, 1, 5, 10} (top to bottom). 



A. Ideal regime 

Consider the ideal operation (CA only) for different 
values of T (assuming the Andreev amplitude to be equal 
to T = 1, and a geometrical factor 7^4 = 0.2). All electron 
pairs are split and cross the device sequentially (one after 
the other). The average (charge) current I (see Fig. 
4a) displays a maximum in the region T ~ -fA T. This 
maximum can be easily explained by the balance between 
the entrance and exit rates of the two electrons in the 
double dot. The total char ge n oise through the entangler 
(see Fig. 5b, and also Ref. Il3j) : 

nch nch 1 nch 1 nch , nch /-i o\ 

D tot — S LL "+" S RR. "T >~>LR ' D RL 

displays a minimum in the same range of parameters, as 
expected. Moreover, for large T/T, it reaches the value 
Ael. This doubling of the Fano factor with respect to the 
Poisson value signals the passage of electron pairs 2 ^. On 
the other hand, the crossed charge noise correlation S L h R 



is positive (Fig. in contrast to the partition noise 

of single electrons, which is negative^. This reflects the 
perfect charge correlation of split pairs 2 ^. As a function 
of T/T, S c L h R also displays a minimum (Fig. (2Ji) and be- 
comes Poissonian in the limit T 3> 7a T, where successive 
pairs are collected much faster than emitted, so are well- 
separated in time. Notice that in the ideal operation one 
has: 

S? R = Sfl = S$ R . (14) 

Turning to the spin noise, one finds that the total spin 
noise S^ t is zero, as discussed above. Moreover, one 
verifies that: 

$ll — S R V R — ~Slr = el , (15) 

thus diagnosing perfect spin correlations from the entan- 
gler. 

B. General regime 

Parasitic processes are next illustrated in presence of 
cotunneling alone (U is chosen as infinite), which intro- 
duces a coupling between branches L, R. In this case, 
pairs are distributed (LL, RR and Li?) but still sequen- 
tially (they pass one after the other and do not mix) : 
therefore all the /'s in Eq. (1) are equal to zero. First, 
assuming Tl = Tr, the average current as well as Sf^ t 
which sums over all branches are not modified 25 . On the 
contrary, S L h R decreases (Fig. and can even become 
negative when cotunneling has a strong effect : due to 
branch coupling, charge cross-correlations are lost and 
partition noise may dominate over pair correlations. On 
the other hand, spin noise gives: 

Str = S s L p R = eIp s , (16) 

Therefore the spin Fano factor simply counts the split 
pairs (Fig. [2)) and measures the entangler efficiency 
One checks that the total spin noise satisfies: 

5^ + 5^ + 25^ = 0, (17) 

i. e. the whole entangler is a perfect Andreev source : 
whatever the origin of the pairs (LL, RR or LR), they 
are well-separated and the spin correlations of successive 
electrons is perfect. Notice that p s tends to one for large 
T S> jcT where CT has no time to occur and becomes 
irrelevant. In practice, for large \E\ — Ei\/T > 5, CT is 
merely suppressed. In summary, with CT alone we have a 
situation where the efficiency of the entangler is reduced 
(and measured by the spin noise), while its fidelity re- 
mains perfect. Yet, notice that further use of these pairs 
in the output quantum channels would require an addi- 
tional device to collect split pairs LR and filter out the 
others. 

Next, double occupancy in the dots is allowed (process 
DA). Fixing T/T, the crossed charge noise is reduced and 
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can become negative at small U (Fig. Spin noise cor- 
relations also decrease as U/T decreases (Fig. |2Ji). The 
structures observed at small U for small T are due to 
the difference between the single electron levels E\, E2, 
and DA transitions becoming resonant when Ei ~ —U. 
The width of the minimum increases for large T, due to 
two competing processes. On one hand, the probability 
for having a split pair (CA) in the 2 dots oscillates on a 
large time scale (•jaT)^ 1 . On the other hand, the proba- 
bility for a DA processes is smaller, but it oscillates more 
rapidly^ ((U 2 + AT 2 )- 1 / 2 ). A high detection rate thus 
dynamically favors the fast DA process, even though it 
is weaker. In the resonant regime Ei ~ —U, Sf^ can 
even approach 4e/L, e.g. noise doubling in branch L. 
Moreover, it is important to point out that due to di- 
rect Andreev processes, the total spin noise is no more 
equal to zero. Actuelly, the occurrence of three-electron 
states in the double dot has the effect of mixing the pairs, 
reducing the fidelity, especially for small U and T. 

In summary, the analysis of this quite general case 
demonstrates that optimal pair correlations are obtained 
if jaT < r < U. Nevertheless, at fixed T/U, a too small 
7A? 1 /r ratio is detrimental, as it favours LL and RR 
pairs by a dynamical effect. 



V. BELL INEQUALITY TEST 

In quantum optics, entanglement is typically probed by 
a Bell inequality (BI) test. It relies on the measurement^ 
of number correlators. One may simplify the response of 
the electronic circuit such as to measure, not the instan- 
taneous current, but instead the particle number accu- 
mulated during a time r is: 



N a (t,r) 



l+T 



dt'l a (t') 



(18) 



(a = ±a, ±6 see Fig. [y. The inequality which is derived 
assuming a product density matrix weighted by local hid- 
den variable readsii: 

|G(a, b) - G(a, b') + G(a', b) + G(a', b')| < 2, (19a) 
((N a (r) - N. a (r))(N b (r) - M 6 (t))> 



G(a,b) 



-(19b) 



((JV (r) + N- a (r))(N a (T) + 2V_ (r))) 

In Refill! the correlators (N a (r)Np(T)) were related to 
zero-frequency current noise correlators via an approxi- 
mate relation. This approximation breaks down for short 
times. We follow Refl2^ and calculate the correlators 
from first principles. The cross correlator for arbitrary 
spin directions of the filters can be separated into a par- 
allel and an antiparallel component: 

(N a (T)N p (T)) = (N L i(t)N ri (t)) sin 2 (0 Q/3 /2) 

+ (A lt (t)A^ t (t)) cos 2 (0 q/3 /2) (20) 

with 9 a p the relative angle between the polarizations in 
L and R. With the same choice of angles as considered 



in Ref. [llj the BI becomes: 

\AN s L y(ANt h R + A+)\<l/V2 (21) 

with 

AA^, = (N L1 (r)N RL (r)) - (N Lr (r)N m (r)) , (22) 
and 

AJV&+A+ = (N L ^T)N Rl (T)) + (N L] (T)N m (r)) , (23) 
where 

A+ = r 2 (I L )(I R ) (24) 
is the reducible part of the charge correlator. 




20 30 



FIG. 3: The left part of the BI is compared to the limit 2" 1/2 
(one assumes T = 1 and 7^ = 0.2). a) It is plotted for 

< Fl.r < 8, two values of E (0, and 5 which simulates 
the ideal working regime), and rT/h = {0.2, 0.5, 1, 5, 10, 20} 
(from the top to the bottom). To characterize the BI test, 
each contribution of | AiV^|/(AjV£^ + A + ) is plotted for r = 
0.2 b), t = 5 c). d) Double occupancy is taken into account, 
assuming E — 5, Fl,r = {0.01, 1,5, 10, 100} (top to bottom) 
and r = 0.2 

In the ideal case (see Fig. 3a), the BI violation is maxi- 
mal for t < h/T (h is the Planck constant). Increasing r, 
thus the number of measured pairs, quantum spin corre- 
lations show a minimum at T ~ ^aT, due to a decrease 
of AN^ ch I A + (Fig. 3b, c). BI violation is recovered 
provided AN^ ~ tI dominates over A + = t 2 I 2 , where 

1 ru j 2 A T 2 /T for large r and I - T for small T. This 
means that pairs should be measured roughly one by 
oneii. Now, considering parasitic processes, CT tends 
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to prevent BI violation at small T, and for large enough 
r (Fig. 3a). Finite U sensibly decreases the quantum 
spin correlations (Fig. 3d). The minimum flattens as Y 
increases, as in Fig. 2d (dynamical effect). Yet, com- 
paring Figs. 2 and 3, one sees that the BI test is much 
less affected by parasitic processes than the spin noise. 
For instance, with T = 5 and U/T = 40, the spin noise 
is about 0.5, signalling a low efficiency, while the BI is 
violated with r = 0.2. In fact, this window allows fil- 
tering of split CA pairs, mostly dropping "wrong" DA 
pairs. Therefore, tuning r to an optimum value yields 
high fidelity of entanglement, even if the efficiency of the 
cntangler is low. 

VI. SPIN-FLIP PROCESSES 

Finally, let us address the effect of spin-flip scattering 
in the dots, here simply described by a spin relaxation 
time r s f. As entanglement is concerned, this acts as a 
decoherence source, as it induces some mixing of singlet 
states (ll s ) with triplet states (11*). Spin-flip is easily 
incorporated into the general density matrix equations^, 
in a way similar to Ref. It consists in defining oc- 
cupation states (n.i CT , n2o-), and adding into the density 
matrix equations the terms corresponding to the decay 
of the spin densities rii a — rii— a . The effect of spin re- 
laxation is summarized on Fig. 4. First, surprisingly 
enough, the average current is strongly quenched as soon 
as the spin-flip rate is of the order of the CA amplitude 
(Fig. 4a). This can be explained as follows. CA is co- 
herent provided it involves a spin-conserving transition 
between degenerate states. Spin-flip in the dot, if faster 
than the CA resonance between states (00) and (ll s ), 
suppresses this resonance, opening a new decay channel. 
Indeed, spin flip competes with CA in a way similar to 
charge decay toward L, R. For t s j < h/Y < Ii/jaT, the 
current behaves as tj*. For the same reasons, spin-flip 
also decreases charge correlations (Fig. 4b) . Secondly, as 
expected, spin-flip decreases the spin noise correlations 
in leads L,R (Fig. 4b), by essentially diminishing the 
fidelity 1 — fg. In addition, Bell correlations are also af- 
fected by spin-flip, but, as above, much less than spin 
noise. 

VII. DISCUSSION 

Let us now discuss the relevance of the above analysis 
in view of a future experimental realization of this entan- 
gler. First, we summarize the required operation regime, 
from the above as well as previous discussionii^ : 

A s , U, |£i - E 2 \ >8e> n LtR , k B &, t~/ 

>Tl,r>1aT, 1c T (25) 

where [Il.r is the voltage drop between the supercon- 
ductor and the leads L, R and Se the separation between 
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r/T r/T 

FIG. 4: (a) Average current I ; (b) spin current noise (con- 
tinuous) and Bell correlation (dotted), as a function of r/T, 
with U — oo,E — 5 (CA only) and finite spin-flip time r s / 
in the dots. r s fT/h — {10, 1, 0.1, 0.01} decreases from top to 
bottom. 

different orbital states in the dots. We assume that the 
average currents in both leads can be measured, as well as 
zero-frequency (charge and spin) noise correlations. Nio- 
bium will be taken for the superconductor (As ~ 9.2K), 
and 6e ~ lmeV ~ 12K for GaAs or InAs-based small 
quantum dots. A reasonable value for U is 45AT, larger 
than As, which does not change the above conclusions, 
for which U < As is instead assumed. A typical relax- 
ation rate is Yl^r ~ lOOmK, larger than ^9. One can 
then discuss the operation of the entangler as a func- 
tion of the CA amplitude, say the parameters T and 7^4. 
As an experimental constraint, we fix the value of the 
current to a minimum of lpA, to make current correla- 
tions measurable. The current is plotted in Fig. 5a in 
nanoAmperes as a function of the electronic amplitude 
T, expressed in Kelvins, for two values of ja- 

T (K) T (K) 



1E-3 0.01 0.1 1 0.01 0.1 1 




0.01 0.1 1 0.01 0.1 



T (K) T (K) 

FIG. 5: Various diagnosis of the entangler, plotted as function 
of the pair tunnel amplitude T (in Kelvins), for ja = 0.01 and 
0.2. a) The average current, assumed larger than 1 pA; b) The 
total charge noise; c) The total spin noise. The spin-flip time 
decreases from bottom to top, T a f(ns) — {100, 10, 1}; d) The 
crossed spin noise. The spin-flip time decreases from top to 
bottom, T s f(ns) — {100, 10, 1}, causing nearly no variation in 
panels a,b. 

For optimum operation, the entangler should provide 
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pairs which are i) well-separated in time; ii) split as much 
as possible in (L,R); iii) in the singlet state. As a first 
probe, the total charge noise tests the temporal separa- 
tion of successive pairs, signaled by the doubling of shot 
noise, S^ t — Ael (Poisson result). We have numerically 
verified (see Fig. 5b) that this criterion is fulfilled up to 
a few percent if jaT < Tl,r. It is compatible with our 
constraint on I for T ~ 0.1K { lA = 0.01) and T ~ 0.01K 
(7A = 0.2). 

As a second probe, the total spin noise also reflects the 
"mixing" of pairs in time, being zero if pairs are well- 
separated (T < r whatever ja without spin flip). Fig. 
5c shows that this criterion is very well fulfilled (lower 
curve). In addition, S^ t is sensitive to spin-flip. There- 
fore a combination of total charge and spin noise gives 
information on pair mixing and spin flip. Spin flip tends 
to increase S^ t from the ideal value 0, especially for large 

1A- 

As a third criterion, crossed spin correlations provide 
information on the fraction ps of split pairs (L,R). We 
find that ja = 0.2 gives an excellent ratio (nearly 1 for 
T < r^j?). On the contrary, for 7^4 = 0.01, this ratio 
drops to about 0.2, LL and RR pairs being "dynami- 
cally" favored. Unless some further processing of pairs is 
achieved, this is detrimental to the entanglcr operation. 

To summarize this analysis, a small value of 7,4 is 
preferable for pair separation, but a large one is re- 
quested for the singlet fidelity. Given these conflicting 
requirements, a good compromise is obtained with large 
ja ~ 0.2 but small T ~ 0.01K. Concerning the spin-flip 
time, values of the order of Wh/T ~ 50ns are suitable. 

VIII. CONCLUSION 

We have analyzed in detail the operation of a realis- 
tic superconductor-dot entangler. We have shown that 



zero frequency charge and spin current correlations al- 
low a detailed analysis of the efficiency and fidelity of the 
entangler in terms of parasitic and spin flip processes. 
Spin Fano factors for spin current correlations appear as 
an optimum characterization of the entangler operation. 
Due to the absence of spin correlation of electrons emit- 
ted within different Andreev processes, such spin Fano 
factors directly probe the spin correlations within an en- 
tangled pair. Depending on the quantity under interest, 
they are equal to ±1 in the ideal operation regime. In 
this sense, measuring the spin current noise is an alter- 
native to time-resolved measurements which should cap- 
ture pairs one by one. Yet, for an absolute diagnosis of 
entanglement, a time-resolved Bell inequality measure- 
ment is necessary, within a narrow enough time window. 
Contrarily to low frequency measurements, this allows 
to check the entangler fidelity even in presence of a low 
efficiency. As a result of this work, we conclude that 
a satisfactory operation is upon reach with realistic pa- 
rameters. For detection rates T ~ lOOmii', the short- 
est time window described for a Bell test corresponds to 
r ~ Ins, a time scale which could be accessible with 
fast electronics apparatus. Moreover, the geometric pa- 
rameter ja appears to be crucial and should be rather 
large. Encouraging results have recently been obtained 
in two experiments which have observed crossed Andreev 
reflections at distances not very small compared to the 
superconductor coherence length^. 
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